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This  report  deal.:;  with  the  current  flow  of  Birkjand 
currents  injected  into  the  auroral  zone  by  solar  events 
and  penetrating  through  the  ionosphere,  stralo-  and 
troposphere  to  the  earth  surface.  In  the  calculat  ion  of 
the  aurora  currents,  the  following  simplifying  assumption:.; 
have  been  made.  The  earth’s  magnetic  field  is  constant 
and  vertical  with  respect  to  the  earth  surface.  The 
problem  is  two  dimensional.  The  potential  function  and 
the  electric:  field  are  functions  of  the  vertical  coordin¬ 
ate  v,  arid  one  horizontal,  coordinate  x  of  a  Cartesian  co¬ 
ordinate  system  x,  y,  z.  The  Hall  current  flows,  only  in 
the  y  direction.  Steady  state  conditions  are  assumed, 
that  is,  the  electric  functions  are  independent  of  time. 
The  conductivity  is  scalar  and  exponential  in  ttie  tropo- 
ancl  stratosphere.  In  the  ionosphere,  the  parallel  conduc¬ 
tivity  continues  to  increase  exponentially,  whereas  the 
Pedersen  or  transverse  conductivity  decrease.';.  These 


conductivity  profiles  are  expressed  by  suitable  analytic 
functions . 

The  earth  is  assumed  to  be  a  plain  equipotontial  sur¬ 
face.  The  current  is  injected  into  the  ionosphere  verti¬ 
cally  in  a  strip  of  7 fi  km  width  and  ejected  in  a  parallel 
strip  of  the  same  width.  The  solution  of  the  differential 
equation  is  given  in  harmonics  or  eigenfunction:-. 


This  solution  may  he  co’isidci  ed  ns,  a  first-  example  to 
show  a  way  how  Lhes"  typo  of  problems  can  ho  solved  anal  y 
licnlly.  The  rest  r  i  c  i  ions  may  bo  relaxed  by  accommoda  t  j  n 
different  conductivity  profiles,  extension  to  lime  depend 
ent  problems,  and  maybe  to  different  coordinate  systems. 

II.  SOL  U  T I  OH  OF  DIFFER  ENT  I  AT,  _F.  QUA  T  ION  . 

The  following  symbols  and  definitions  are  used  in  thi 
calculation : 

x,  y,  z  =  Cartesian  coordinate  system  with  the  z  axi.s 
pointing  upwards. 

i,  j,  k  =  Unit  vectors  in  direction  of  the  coordinates 
>>1  -  Xg  exp(2kz)  -  scalar  and  parallel  conductivity 
Aq  =  2  x  10  ^l/Uni  =  conductivity  at  the  ground 

2k  =  !”  --0  -  scale  factor  of  conductivity 
“  1 0  k  m 

2 

X ^  ~  A l/ ( 1 i a  exp (2k z)  )  =  Pedersen  or  transverse 

conduc t i v i ty . 

t t>  =  <}>(x,z)  -  Potential  function 
E  --  —  =  electric:  field 

T  =  L  E  (I,  is  the  conductivity  tensor  in  the  ionosphere 
and  a  scalar  in  the  strain-  and  troposphere) 

Fl  =  magnetic  field  vector 
)j  -  permeability. 
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Since  thi-  current  flow  .it;  assunic'!  to  bo  in  i  early  r;l  at 
from  (2) 


-  0;  V  X  K  0;  ]■:  -  v;  (1) 

ot 

and  taking  the  cl  ivory  once  of  (1) 

V  •  V  x  fi  =  0;  V  •  T  ~  -  \-  (I\2  i?:  -I-  :jA  ?(l4^-x  13 )  /IS 2 
+  -  <>  <4> 


The  divergence  of  the  Hall  current  in  (4)  is  zero  because 
it  requires  differentiation  to  y  and  the  Hall  current  of 
our  problem  is  not  a  function  of  y.  Executing  the  diver¬ 
gence  operation  in  (4)  we  obtain  the  differential  equation 
of  our  problem 


c)  ,  ou 

iix  2  ';>•< 


a__ 


c'Z 


0 


Since  A ^  is  a  function  of 
d  if  C e r  e n  t  i a  1  ope r  a t or  to  x 


only, 
a  i  i  cl  d  j 


(5) 

we  may  put  .it  before  the 
ide  ( 1)  by 


This  form  of  the  d  i  1' !  er  <.  nt.  j  a  1  equal  ion  <  (,) 


i  mmed  ia  tel  y  a  solution  in  the  form  of  e  inenf  imet  ions  mi 
harmonics.  We  assume:  that  •;  j given  a:;  a  jm  o  !u«-l  of  a 
function  f  =  f(z)  and  q  -  q  ( ) 


<l>  ~  g  ( x )  f  ( 7. ) 


By  writing  the  factor  not  affected  !>y  t!:e  d  i  f  (ercnl  i  ation 
before  the  differential  operator  and  dividing  by  <;■  splits 
(6)  into  two  terms,  the  first  depending  only  on  x  and 
the  second  on  7. 


I  <L_a  +  ...i_  <L.X  --  .o  (8) 

g  dx'  A2?  dz  dz 

We  set  now  the  first  term  equal  to  the  arbitrary  constant 


1  d2q 


with  Lhe  solution 

g  =  A  oxpflisx)  ;  i  ~  /-.I  (10) 

A  is  here  a  complex  integration  constant  and  s  a  scaling 
factor  of  x. 

Replacing  the  first  term  in  ( 8 )  by  --s'  we  obtain  foi 
f  the  differential  equation 


d  .  d  f  2 

-.--A.  -  -  s  A  f  -  0 

d  z  1  d  z  2. 
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w-xp  2ky.) 


f  ■-  0 


Th°  fum;tion  clropj;  out  .uid  wo  r:viy  i  f  v  it 

/  ^  j  in[j.L  ]  l  y  t  ho  equation 


furthor  hi'  the-  auhstituiion 


exp(~2kx); 


rtu 


dz 


--  ,  k  Li 


( 12) 


Wo  obtain  then  t.ho  simplc.  rtj  i  ror,.ft,  iCi 
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c-qua  tion 
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b’o  sot:  now 


f  :=  ( u  i- a ) 


.. ,  m 


(14) 


rt2f 
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(15) 


anrt  «u!«;tituto  c->:jh  ckj  i  on  (10)  in  (U).  This  gi 
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III.  DISCUSSION  OF  WK  SOLUTION 

The  functions  <;ij  and  $  represent  the  in  and  out  nolo 
functions  of  the  differential  equation,  meaning  that  the 
inpole  function  ^  is  caused  by  a  source  at  the  ground 
(•/.  -  0),  where  1  ^  has  its  maximum  value  that  decreases 
with  increasing  distance  from  the  source.  The  outpole 
function  is  produced  by  a  source  at  some  altitude  H 


f> 


. i i  ' > v  •.ant! ,  for  j  nsi  ai  'he  i oic  mpher or  .,1  >v  , 

has  _j  •  ;  r  i  i ;  its  max  i  .!.  J.T1  Vill  !••-•  Lit  th  •  ■  (7,  -  li)  that 

deer  -a  i  I  a  i  1  n :  1  im  s  i  ncj  di.-;  .nice  1.1:0.;:  the  Soo'er,  i  . . ;  . 

a:;  wo  ..pp  roach  the  ground  (a  0).  Thi.;  d  is  t  i  n  •  t  ion 
!'(■!  in  and  out  polo  potent  in.1  fund  ions  is  reflected 

in  modern  terminology  an  up  (.source  below)  or  down  ( source 
above)  inapp  ing . 

To  e.-cc.l  udo  the  ionospheric  and  restrict  t  he  solut  ion 
to  the  tropo-  and  stratosphere  where  A  j  -  A.  wo  rot  a  0. 
The  solution  for  cj  is  not  affected  and  remains  the  same  as 
given  in  (10).  However,  fj  and  f7  reduce  to  the  simpler 
f  o  rm 

-in. 2kz 

f1  -  e  1  (23) 

-no  2k  2 

f?  -  e  z  (24) 


resulting  in 


-  A  oxp('isx)  exp(-(v  1  1  (r,/k)  '  +  1)1:2) 


/;■ 


;>:)  exp  (  (VI  (s/k 


/k)  - 1 )  k '/ ) 


(20) 

(20) 


1(1  p  --  A  exp  (  t  is: 

In  this  reduced  form  it  is  easy  to  recognize  that 
cr‘>  ^  is  the  inpolo  function  with  its  maximum  value  at  the 
ground  z  -■  0  and  diminishing  with  increasing  z .  s  ■  2  f; 
the  outpolc  function,  that  has  its  minimum  value  at  the 
ground  and  increases  with  increasing  /.  Furthermore,  we 
see  that  the  inpole  potential  decreases  much  faster  than 
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t.iic  on  tonic  pot  on  1  ul  Lucre.  ::e:;  .  >.  !•’  J.<  >  i  [•  )•,  < . 

Ihat  t  he  a  t tenant  ion  bv  »n j  .i  ng  down  if.  rush  smu.i.  !  1:  than 
tlie  ci  I  V euuat  ion  b-r  mapping  up,  ;i  f  .  •:/ f  '..hie!,  in  often  u. io 
od  .in  the  litoralure.  The  can.  i.y  t.  i  r  - ,  s  |  solution  <\  Ives  the 
additional  information  t  ha  I  the  v;oak  attenuation  by  down 
mapping  depends  also  on  the  ratio  of  s/k.  Ij  thin  ratio 
is  much  larger  than  1  meaning  for  vide  spread  source.-,  the 
attenuation  by  up  and  down  mapping  will  be  of  the  same  aide 
We  may  reduce  (2b)  and  (26)  further  by  multiplying  the 
bracket  in  the  .second  exp-function  with  k  and  then  let 
k->0 .  This  would  reduce  the  conductivity  to  a  constant  and 
and  (>2  to 

-  A  cxp(tisx)  exp  ( —  s /. )  (27) 

('2  =  A  o  x  n  ( -i  i  s  >: )  c  x  p  ( +  s  x )  (2  fj ) 

These  equations  give  the  solution  for  the  electrostatic 
problem,  a  solution  that  wo  will  need  for  solving  time 
dependent  problems.  Here  the  decay  for  up  and  down  mapping 
is  the  same.  This  shows  quite  clearly  that  the  difference 
in  the  attenuation  in  the  previous  cases  is  caused  by  the 
chang i ng  conduct i v i ty . 

IV.  APPLICATION _0F  THU  SOLUTION  TO  AN  _10N0dPHKR  TC  _PHOHI.EM 
l’or  the  application  of  eigenfunctions  of  a  differential 
equation  to  a  specific  physical  problem,  ww  have  to  specify 
th'>  liouadary  condition  and  l  he  source  of  flit'  potential 
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funo!  .ion  .  I  Km  o  v/»  *  th--'  !•’!  :  ow ■ ng  bound. ;i  y  c  aid  ■  L  ;  *ou 

I.  'i'li  u  cur  i  ont  f  luv;  i  : .  con!  in.  •  !  l<,  a  n  in!  i:\iic  ]  onq  cvl  i  n- 

d<  r  \:  i  t'  h  a  t ; i :  1 1 : ,  .  > !  i  ■  r ,  ■  .  soot  ion  or.  .''how;!  in  ]•'  i  g .  I  . 

..  _  r  ') 

■  Cun  with  an  ivcniy  .•nrrniL  rh.os.i  :.y  of  .1  ].(j  /,/<:" 

is  injected  into  tin  1.  cl  1  hull  of  th-'>  oyi  indor.  Lop  no! 
eject  ed  out  of  llin  right  hull  of  the  top  .air  face. 

3.  No  current  should  penetrate  the  :  >.do  walls  of  the  cyl  i n 
dor.  For  x  r-  0  and  x  150  km  «.,'••/ =  0. 

4.  The  bottom  plain;  (earth  uur  face)  shall  be  an  to  poten¬ 
tial  surface  with  the  potential  value  0. 

5.  The  parallel  conductivity  and  the  Pederson  c  '  vtiv 
ity  A  r)  are  given  by 


*]  *  >  ()  exp  ( 2kz) 

2 

A  2  =  A  i  /  ( .1  +a  exp  ( 2  k  x )  ) 

The  solution  of  the  differential  equation 
state  conditions  has  been  given  in  the  font  of 
tions  in  section  IT,  equations  (21)  and  (22). 


(4.1) 


for  steady 


o  ig  on  f  'me¬ 


lt  remains 


to  construct  by  superposition  of  these  oigonf miet .i ons  with 
suitable  adjusted  constants  A  a  solution  that,  fits  the 
boundary  cond i t ions . 

First  we  accommoda te  the  boundary  condition  3.  for  the 
inpolo  function.  ft/  cix  -=  0  for  x  -  0  and  x  =  150  km. 

From  (21) 


kf  j  /  'A;  ~  A  is  exp(tisx)  f  ^ 


Asf^  i (cos ( sx) 


i  sin ( :  x) 


i 
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c  ;  .  \  r  \ 


!>  I.  Of  {  '  .  ')  ?  '  • 

I  I  ■. 

i  •  0  i'vi'l  ' 


'i'li -  :-r c  fore 


f.  s in  ( 


(  A  .  • ) 


and 


-  A  f1  cos 


(A  .  5) 


The  same  boundary  condition  applies"  also  In  the  out  polo 
function,  r  osu.l  t  iwj  in 

v,  -  I!  r,  cos(sx)  (4.6) 


w.i  tli 


c,b  •  . 

sin's:;)  (4.7) 


To  on  Tore 

o  the?  boundnr 

y 

•ond 

it  ion  4  ,  <; 

0 

for  z  =- 

0  , 

apply  1  ho 

it1  i  r  r  or 

1  nic’ujo 

n;e  l 

hod 

and  sup'-ri 

i:ino.’;o 

«;•  1  and 

'*  2 

*  =  •'!  +  ■‘Z 

-  (A  r 

s  :  11 

r.,) 

i. 

co. 

s.  (s,  ) 

(4 

0 

for  x  -  0  follows  iron  (1.9)  and  (20) 

—  -  “  '  l  » 

,}  (  /  I  (  s/k )  2  VI) 

f ] 0  •  n  + 

.  J  (  /l  I  (s/k)  2  -  J  ) 


I 


I 


2.0 
I  0 


(4.10) 


and 


^  R^20  ^  ^  20 


f  l/f  1q  )  c°;j  (s>: )  -  B  F(z)  cos  (;•.:•:)  (4  .  1  1) 


wi  th 


F  ( z ) 


f 


20 

10 


(4.12) 


The  last  boundary  condition  to  bo  fulfilled  is  condition  2. 
To  do  this  wo  have  to  calculate  the  average  current  density 
at  the  source  for  z  ~  H. 

We  set 


M 

z= II 


!•”  (H) 


1(H)  =  i  A  2Hs  F'xin  sx  +  k  }  B  F'(F1)  cos  r.x  (4.11) 

Since  at  the  altitude  of  z  =  II  -  150  krn  the  Pedersen  conduc¬ 
tivity  A ^  is  very  much  smaller  than  the  parallel  conductivity 
X j  we  may  neglect  the  horizontal  against  the  vortical  current 
flovr  component  and  write 

I  (II)  -  A  j  B  F '  (II)  coy.  (r.x)  (4.14) 

This  approximation  '  is  an  expression  of  the  magnetic  field 
aligned  current.  The  magnetic  field  is  vertical  and  there¬ 
fore  only  the  vertical  current  component  is  significant. 
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The  constant  R  i.e  in':  .  i  c'  r  i:;  i.  rod . 

n  ,  *  bl  ( 

A-,  I"  (!i)  ?  •  {  !■•  (II) 

V  .  Nl  I  HR  !<  1  CAL  KVA  i  .1  IAT !  OK 
.1 .  Conduct  ivitioj:  A  (  anti  ' 

A-j  ~Afi  exp  (2k/) 

>  =  2  X  10“ 11  1  /■'!.! 

?k  .  In  10 
1.0  km 
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run  id 

t  1 1 0 
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-  20  1  /f ’in - 

A .  " 

A 

*0 

exp.  (.2 

( 1.  ha 

k  / ) . 
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(  h  .  1  ) 
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Ai ’oroiicl .  i  ivj  t  ho  lower  tonospl,  or  o  A 


.1  :n 


alee  i  km: 


'  1 


v; i  t  h  in.rvc  i r,  in..:  !  '  i  :  a 


iiniim  .siui  f  i  ns !  1  y  decays 
TIik-  max  i  mum  if.  reached  when 


a  -  exp  (-21:7.  ) 


(5.2) 


If  wo  would  like  the  maximum  of  A ^  to  he  at  100  km 
altitude  then 
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Fig.  2  shows  the  conductivities  of  A^  and  A^ 
given  on  the  horizontal  axis  with  a  logarithmic  scale 
as  a  function  of  altitude  z  on  the  vertical  axis.  The 


maximum  of  i-s  a,:  altitude  with  a  =  /  10  x  10 


A 


-8 


2.  Potential  Fund  ion 


(|>  =  B  F(>i)  cos(s.v) 

B  =  -  if  I  /  2  A  j  ( H)  F '  (11)  -  11.4  MV 

Fig.  3  shows  a  plot  of  the  oqu  ipoten  t  i  a  .1  lines,  with  the 
zero  potential  as  a  straight  vertical  line  in  the  middle 
connected  to  the  ground.  The  lines  to  the  left  increase  in 
steps  of  about  1  MV  and  the  lines  to  the  right  decrease  in 
steps  of  about  -1  MV.  It  is  remarkable  that  the  equipoten- 
tial  lines  in  the  ionosphere  from  1130  km  to  about  75  km 
altitude  arc  practically  vertical.  Then  they  slant  slightly 
to  the  left  in  the  left  half:  and  to  the  right  in  the  right 
half  of  Fig.  3  until  they  bend  sharply  in  a  horizontal 
direction  as  they  approach  the  ground.  This  leads  to  a 
crowding  of  the  equ  i  poten  t  i  a  1  lines,  in  the  surface  layi  of 
10  to  20  km  thickness  and  oilers  the  opportunity  to  measure 

1.3 


the  p'd  end  i 1  dial'  i  !>■•:'  ion  in  ‘  !;<  i  oni  i  bo  i  • i 

a  i  r;d  ;mf'  sounding  of  t  h  :  a  l inocii1!!  t  i.<:  el  <:  i  c  ;.u.  on*.  ini", 

in  i  hr*  lowest  layer  of  !  h*  a !  .fosnh'M  o . 

Func!  ioi;  <!>/•;,;  v^rncr;  a]  t.  i  '•  mi  •  7  in  plot  l  od  in  Fig. 4. 
Tho  horizontal  axis  r;li  §/  from  0  Lo  i  and  t  ] : •  •  \  ei  t  ica 
axis  the  altitude  /  froni  0  to  InU  km.  In  10  kr.’.  a  It itude 
we  have  already  reached  00  5  of  «-he  final  value  ol  V.  The 
shape  of  the  curve  is  the  same  for  any  distance  from  the 
midpoint  x  -  7  5  km,  only  the  f  inal  value  changes  with  co«  (s 
from  the  maximum  value  of  about  +10  MV. 


3.  Electric  Field  Components  E  and  E? 

E  -  s  D  F(z)  sin  sx 
E?  -  -B  F’ (z)  cos  sx 

The  magnitude  of  the  electric  field  components 

F.  and  E,  versus  altitude  are  shown  in  Fig.  5.  The  hor .i- 
x 

zonta.l  axis  displays,  the  fields  E  ,  and  the  vertical 
axis  tho  altitude  7.  The  vertical  field  E  has  its  maxi¬ 
mum  at  the  ground  of  about  ‘2.3  kV/m  and  decreases  rapidly 
with  increasing  altitude.  At  10  km  altitude  the  value  is 
E  =  14  0  V/ni. 

The  horizontal  field  E  is  zero  at  tlie  ground  reaches 
200  V/rn  at  10  km  and  about  220  V/m  at  20  km  altitude.  From 
here  on  it  keeps  its  value  all.  the  way  up  to  150  km  altitude. 
This  is  indeed  a  very  remarkable  result..  Since  the 
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Another  dishing  uinh  iiKj  character  ist  ic  in  the  phase 
shift,  between  the  horizontal  and  vertical  component  K  and 
E^ .  For  x  =  0  and  x  -  ISO  an.  the  ver  t  ical  component  E 
has  its  maximum  or  minimum  value  of  '750  V/m  at  10  km 
altitude.  At  x  -  15  km  it  passes  through  zero  and  changes 
polarity.  The  horizontal  component.  E  is  zero  at  x  =  0 
and  x  =  ISO  km  and  has  its  maximum  value  of  200  V/m  at 
x  =  75  km  where  Ev  passes  through  z^ro .  An  airplane  flying 
at  10  km  altitude  from  x  =  0  to  -  150  km  would  record  a 
field  pattern  as  shown  in  Fig.  6.  The  horizontal  axis  gives 
here  the  flight  path  x  at  10  km  altitude  and  the  vertical 
axis  the  field  components  E  and  E  .  E  is  displayed  by  a 
dashed  and  E  by  a  solid  line.  It  is  most  unlikely  that  such 
a  field  pattern  could  be  caused  by  an  atmospheric  electric 
fair  weather  field.  The  constant  fair  weather  field  would 
appear  as  a  zero  shift  of  the  Ez  component  which  could  easily 
be  eliminated  from  tfyo  record. 
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Fr  on 

75  Lo  100  km  altitude  the  horizontal  current  density 
decreases  again  and  becomes  ncyl.  iyable  in  altitudes  of  100 
to  150  km.  The  vertical  eomnonent  1  starts  to  increase 

‘  7. 

from  about  55  km  reaches  approximately  its  maximum  value 
-  5  2 

of  1.37.10  A/m'  at  100  inn  and  stays  at  this  value  up 

to  150  km  altitude  (field  aliened  current). 

Fiy.  0  gives  the  cunenl  profiles  of  Ix  and  I  in  the 

lower  layer  of  the  atmosphere  from  0  to  30  km  altitude. 

- 1 1  ? 

I  is  4.7x10  A/in  '  at  the  ground.  This  is  already  higher 
by  about  a  factor  30  than  the  fair  weather  current.  It  has 
doubled  its  value  at  about  20  km  altitude  and  increases 
more  rapidly  with  further  increase  of  altitude.  The  horizon¬ 
tal  component  I  starts  with  zero  at  the  ground  reaches  the 
value  of  I  at  about  11  km  altitude  and  increases  much  more 
rapidly  than  I  for  higher  altitudes.  Note,  however,  that 
that  phase  shift  between  Ey  and  I-:  discussed  above  applies 
also  to  I  and  I  r/ . 
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It  should  not  I .  ■  difficult  to  design  front  the  Ch  irac!  *•!.- 
istic:  dif f eronco.;  in  the  field  .mk!  currcnl  jMtlerii  oi  fin- 
fair  weathei  and  aurora  current  flow  an  evaluation  method 
f ot  the  data  obtained  from  a  network  on  the  ground  or  from 
an  airplane  to  separate  the  fair  weather  from  the  aurora 
fields  or  currents.  It  is  obvious,  however,  from  Pig.  7  and 
Pig.  8  that  the  current  density  1  of  the  aurora  currents  is 
not  constant  with  altitude.  Since  a  constant,  current  density 
with  altitude  is  a  necessary  assumption  to  extrapolate  the 
fair  weather  potential  at  flight  altitude  to  the  ionosphere 
to  obtain  the  ionosphere  potential  of  the  atmospheric  electric 
current  the  same  extrapolation  method  would  not  be  valid  to 
determine  potential  differences  in  the  ionosphere  potential 
caused  by  solar  or  ionospheric  events. 
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